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Altered gene expression in dry age-related macular degeneration 
suggests early loss of choroidal endothelial cells 
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Research, The University of Iowa, Iowa City, lA; ^Department of Biomedical Engineering, The University of Iowa, Iowa City, lA 

Purpose: Age-related macular degeneration (AMD) is a major cause of blindness in developed countries. The molecular 
pathogenesis of early events in AMD is poorly understood. We investigated differential gene expression in samples of 
human retinal pigment epithelium (RPE) and choroid from early AMD and control maculas with exon-based arrays. 
Methods: Gene expression levels in nine human donor eyes with early AMD and nine control human donor eyes were 
assessed using Affymetrix Human Exon ST 1.0 arrays. Two controls did not pass quality control and were removed. Dif- 
ferentially expressed genes were annotated using the Database for Annotation, Visualization and Integrated Discovery 
(DAVID), and gene set enrichment analysis (GSEA) was performed on RPE-specific and endothelium-associated gene 
sets. The complement factor H (CFH) genotype was also assessed, and differential expression was analyzed regarding 
high AMD risk (YH/HH) and low AMD risk (YY) genotypes. 

Results: Seventy-five genes were identified as differentially expressed (raw p value <0.01; >50% fold change, mean log^ 
expression level in AMD or control > median of all average gene expression values); however, no genes were significant 
(adj. p value <0.01) after correction for multiple hypothesis testing. Of 52 genes with decreased expression in AMD (fold 
change <0.5; raw p value <0.01), 18 genes were identified by DAVID analysis as associated with vision or neurologic 
processes. The GSEA of the RPE-associated and endothelium-associated genes revealed a significant decrease in genes 
typically expressed by endothelial cells in the early AMD group compared to controls, consistent with previous histologic 
and proteomic studies. Analysis of the genotype indicated decreased expression of ADAMTS9 in eyes with high-risk 
genotypes (fold change = -2.61; raw p value=0.0008). 

Conclusions: GSEA results suggest that RPE transcripts are preserved or elevated in early AMD, concomitant with loss 
of endothelial cell marker expression. These results are consistent with the notion that choroidal endothelial cell dropout 
or dedifferentiation occurs early in the pathogenesis of AMD. 



Age-related macular degeneration (AMD) is the leading 
cause of blindness among the elderly in developed countries 
[1]. AMD involves the progressive loss of photoreceptor cells 
from the macular region of the retina, resulting in impaired 
vision and, in advanced stages, blindness. At least three cell 
layers undergo changes in AMD, including the photoreceptor 
cells, retinal pigment epithelium (RPE), and choriocapillaris. 
The RPE regulates the activities of the photoreceptor cells and 
choriocapillaris. For example, RPE cells actively phagocytose 
photoreceptor cell outer segments, recycle vitamin A, shuttle 
debris from the photoreceptor cells to the bloodstream, and 
import glucose, oxygen, and other components to accommo- 
date the high metabolic demands of the retina [2], in addition 
to providing trophic support to the choriocapillaris [3,4]. The 
choroid serves as a high-volume transportation courier, deliv- 
ering nutrients to the RPE and accepting waste products for 
further processing elsewhere in the body. The preclinical and 
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early stages of AMD are recognizable by increased forma- 
tion of lipid-rich sub-RPE deposits termed drusen and altered 
RPE pigmentation [5,6]. 

The photoreceptor cells, RPE, and choriocapillaris 
endothelial cells form an interdependent complex. Injury or 
dysfunction in any of these layers leads to loss of the other 
two in several chorioretinal diseases. A more complete under- 
standing of the early sequelae of events in AMD is necessary 
to guide new therapies. Numerous interdependent biologic 
processes have been implicated in the pathogenesis of AMD, 
including increased activity of the complement cascade, 
infiltration of cells mediating inflammatory responses, 
increased oxidative stress, and altered lipid metabolism [7,8]. 
Although RPE cells are typically viewed as the primary cells 
affected in AMD, changes in the microvasculature of the 
choroid (choriocapillaris) have also been reported in asso- 
ciation with drusen, including dropout of vessels [9,10] and 
decreased blood flow [11]. In a subset of advanced AMD 
cases, choroidal neovascular membranes (CNVs) form as 
blood vessels from the choroid breach the RPE and proliferate 
either beneath the RPE or in the sub-retinal space. Expression 
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Figure 1. Overview of bioinfor- 
matics pipeline. Software is shown 
in the upper left corner of the boxes 
while the analytic process is indi- 
cated in the center of the boxes. 



of vascular endothelial growth factor (VEGF), a marker of 
hypoxia, has been implicated in the formation of CNVs [12]. 
In current medical practice, only after CNVs have appeared 
and photoreceptor cell death has occurred can therapeutic 
measures be taken to slow further vision loss [13]. 

Despite considerable progress in unraveling genetic risk 
factors for AMD, major challenges remain. The relationships 
between the biologic processes remain uncertain, and the 
initial molecular conditions driving development of AMD are 
poorly understood. Evaluating gene expression in early AMD, 
intermediate AMD, and advanced AMD is one approach to 
advancing exploration of these problems. The first large- 
scale study of gene expression in the AMD-affected retina 
and RPE and choroid tissue identified changes between 
various stages of AMD, including apoptotic and neovascular 
pathways in advanced AMD [14]. As part of a study exam- 
ining the relationship between AMD and gene methylation. 
Hunter and colleagues examined gene expression in AMD 
and normal samples [15]. They found that expression of 
glutathione S-transferase isoform mul (GSTMl) and mu5 
(GSTM5), antioxidant isoenzymes, was reduced in the RPE 



and choroid of eyes affected by AMD compared to control 
eyes; this reduction was correlated with hypermethylation of 
the GSTMl promoter [15]. However, only two of the samples 
were classified as early AMD, and analysis was performed 
without respect to AMD grade. 

To further investigate the molecular events initiating 
AMD, we evaluated gene expression in early AMD and 
normal RPE and choroid tissues using exon-based micro- 
arrays (Affymetrix). Our analysis identified a statistically 
significant decrease in expression of endothelial genes in 
early AMD with preservation of RPE-specific transcripts, 
suggesting that vessel loss or dedifferentiation may precede 
advanced damage in AMD. 

METHODS 

Tissue acquisition: Donor eyes were obtained through the 
Iowa Lions Eye Bank (Iowa City, lA) after receiving informed 
consent in accordance with the tenets of the Declaration of 
Helsinki. Donor age, sex, cause of death, and diagnostic status 
are listed in Table 1. Eyes were dissected, and 6 mm punches 
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Figure 2. Quality control plots generated by arrayQualityMetrics before and after removal of CTRL3 and CTRL9. A and B are a false color 
heatmaps indicating the distances between arrays, computed as the absolute mean distance between the data. C and D indicate the sum of 
distances computed for A and B. For C the outlier threshold was 7.14 (vertical bar), and for D the outlier threshold was 5.79 (vertical bar). 
When analyzing all 18 samples, both CTRL3 and CTRL9 were flagged as exceeding the outlier threshold (A and C). When CTRL3 and 
CTRL9 were omitted, only ARM7 exceeded the threshold (D). E and F are boxplots of the gene-level signal intensity distribution across 
the arrays. For E and F, the Kolmogorov-Smirnov statistic was calculated based on the distanced between each individual array and the 
pooled distribution of all arrays. G and H show for each array with outlier thresholds (vertical bars) of 0.0509 and 0.049, respectively. 
Only CTRL3 was flagged by this metric (G). 



were taken from the macula, centered on the fovea centralis. retina, flash frozen in liquid nitrogen, and stored at -80 °C. 
Macular RPE and choroid were separated from the neural The neural retina was used for separate experiments. In this 
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Figure 3. Distribution of mean 
expression before median thresh- 



2 4 6 8 10 12 olding. Mean expression is log^ 

of gene-Ievel signal intensity. The 
mean expression median indicated by the red line. 



manner, RNA was stabilized within 4 h of death. RNA was 
extracted from the frozen RPE and choroids using a commer- 
cially available kit (RNeasy Mini Kit; Qiagen, Valencia, CA). 

CFH risk-allele genotyping: To characterize how the risk 
allele in complement factor H (CFH; single nucleotide poly- 
morphism [SNP] rsl061170) affects gene expression in the 
RPE and choroid, genotyping was performed on DNA from 
either whole blood (collected in EDTA-coated tubes) or extra- 
ocular muscle. DNA was isolated using established methods 
with the Centra system (Qiagen, Valencia, CA) or the Qiagen 
DNeasy Blood & Tissue Kit, for blood and muscle, respec- 
tively. Cenotyping of this polymorphism was performed 
using TaqMan predesigned SNP genotyping assays (Applied 
Biosystems) in a high-throughput system (Fluidigm, San 
Francisco, CA). 

Microarray processing: Exon expression levels of the 
RPE and choroid RNA were determined using Affymetrix 
GeneChip Human Exon 1.0 ST arrays processed in two 
batches at the University of Iowa DNA Facility. Samples 
were used essentially as previously described [16]. Briefly, 
single primer isothermal amplification was used to convert 
25 ng total RNA to cDNA with the WT-Ovation Pico RNA 
Amplification System (NuGEN Technologies, San Carlos, 
CA). The resulting cDNA was purified using a Qiagen 
QIAquick PCR Purification column, converted to sense target 



(ST)-cDNA using the WT-Ovation Exon Module vl (NuGEN 
Technologies), and purified once more. ST-cDNA was then 
fragmented (85 nt mean length), and the NuGEN FL-Ovation 
cDNA Biotin Module, v2 (NuGEN Technologies) was used 
to biotin-label the fragments according to the manufacturer's 
guidelines. This product was combined with Affymetrix 
eukaryotic hybridization buffer (Affymetrix, Santa Clara, 
CA) and hybridized to Affymetrix Human Exon 1.0 ST arrays 
(Affymetrix). An Affymetrix Model 3000 scanner with 7G 
upgrade was used to scan the arrays. Data were collected 
with GeneChip operating software (ver. 1.4; Affymetrix) and 
saved as CEL files. 

Microarray analysis: 

Preprocessing, expression analysis, and quality 
control — ^To assess differential expression and alternative 
splicing, CEL files were normalized by robust multiarray 
average (RMA) [17] to facilitate comparisons across arrays 
using AltAnalyze (ver. 2.0.7 beta) with gene and transcript 
data from the Ensembl 65 database and human genome 
build Hgl9 [18]. RMA, a common algorithm for processing 
microarray data, corrects for background fluorescence, 
normalizes data for comparison between arrays, and produces 
logj-transformed estimates of probeset expression. Technical 
variation among microarrays is well known and often arises 
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Figure 4. Heatmap of differentially expressed genes shown in Table 2. Dark shading indicates low expression; light shading indicates high 
expression. 



when arrays are processed in separate batches (i.e., prepared 
by different people, on different days, etc.) [19]. To address 
this issue in our data set, initial probeset-level expression esti- 
mates were corrected for batch effects while preserving AMD 
status, CFH genotype, and sex covariates using the ComBat 
algorithm [20]. ComBat removes batch effects more effec- 
tively than several other procedures [21]. The batch-corrected 
data were reimported into AltAnalyze and processed as 
follows: masking of cross-hybridizing probesets; filtering 
probesets expressed below non-log level (i.e., expression 
intensity before log, transformation) of 70; gene-level differ- 
ential expression using a moderated t test; and computation 
of false discovery rate (FDR) adjusted p values to account for 



multiple hypothesis testing correction [22]. Heatmap genera- 
tion and hierarchical clustering were performed using the R 
statistical software (ver. 3.0.0) [23]. To cluster genes based 
on the similarity of the expression pattern across samples, 
we used a Pearson-based distance metric (1 minus Pearson's 
correlation coefficient) [24]. To assess alternative splicing, 
the splicing index and MIDAS [25] procedures in AltAnalyze 
were used. Putatively alternatively spliced transcripts were 
visualized using DomainGraph (ver. 3.0) [18], a plugin for 
Cytoscape (ver. 2.8.1) [26]. Quality control metrics, including 
distance between arrays and comparison of array intensity 
distributions, were calculated using the arrayQualityMetrics 
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package (ver. 3.16.0) [27] for R. After outlying arrays were 
removed, the final AltAnalyze results were recomputed. 

Gene set analysis: The Ensembl BioMart tool was used to 
map identifiers from previous gene expression studies to 
Ensembl IDs [28]. Gene set enrichment analysis (GSEA; ver. 
14) [29] implemented with GenePattern (ver. 3.6.0) [30] was 
used to evaluate the overrepresentation of custom gene sets 
with AMD or control samples. For GSEA, phenotype permu- 
tation was performed using 1,000 permutations and two gene 
sets, an RPE-specific set and an endothelium-associated set, 
compiled based on literature search. Annotation of differen- 
tially expressed genes was performed using the Database for 
Annotation, Visualization and Integrated Discovery (DAVID) 
Functional Annotation Tools with default parameters [31]. 

RESULTS 

Quality control: We defined early AMD as RPE changes 
(depigmentation or hypopigmentation) and/or macular drusen 
without geographic atrophy (GA) or CNV [32] as determined 
by chart review. Eyes with CNV or geographic atrophy were 
excluded. Since early AMD was classified in the donors by 
ophthalmoscopy, completely ruling out a rare event such as 



occult CNV is impossible in the early AMD samples. With 
the exception of one donor who later did not pass quality 
control, all unaffected control donors were required to have 
had a normal fundus exam within two years of death. Based 
on these criteria, nine samples were classified as early AMD 
(designated ARM 1-9) and nine samples were classified as 
controls (designated CTRLl-9; Table 1). To assess the statis- 
tical equivalence of donor ages in each group, we performed 
Welch's two-sample t test, which gave a p value of 0.54. The 
mean RNA integrity numbers for samples processed in batch 
A was 6.76 (standard deviation [SD] 0.34) and for batch B 
was 6.22 (SD 0.59; Table 1). Following normalization in 
AltAnalyze (Figure 1), samples CTRL3 and CTRL9 were 
flagged as potential outliers by arrayQualityMetrics analysis 
[27] based on their sum of distances to other arrays (Figure 
2A,C). Additionally, the gene-level signal intensity distribu- 
tion of CTRL3 was lower than that of other arrays and was 
flagged by arrayQualityMetrics's Kolmogorov-Smirnov- 
based outlier detection module (Figure 2E,G). Given the 
magnitude of the difference between these two arrays and the 
other arrays, we removed these two arrays from the data set 
and reprocessed the remaining 16 samples (Figure 2B,D,F,H). 
With this reduced data set, ARM7 was flagged as an outlier 



2284 



Molecular Vision 2013; 19:2274-2297 <http://www.molvis.org/molvis/vl9/2274> 



) 2013 Molecular Vision 



PS S 



PS 



en 



S < 



S 



o 



s 



O 
(N 



Q 



P s 

T3 I- 



a 00 



o 



o 



o 

1 1 

60 



■S 



O 



60 
O 
J3 
& 

O 

S 

tu 
c 
o 



60 



c 
a 

60 
O 



3 ^ 



60 
C 



S 



Pi 



J3 



tu 



o 



o 

J3 



J3 

n, 

O 



3 tu .ti _^ 



tj 



.S^ " B 



"5 m 



Q 



o 

O 



a 

60 

c 
c 

■a 
■§ 

■& 

o 



o 

D. 
60 
C 

■3 
c 



60 

c 

'-3 
a 



• — ^ „ a J2 



„ 3 

IS tu 
c S 



-^3 



CO 



tu 



tu 



tu 



o 



s -2 .a 



to 





00 


t3N 














ta^ 


1^ 














1^ 


<ri 


m 




-a- 




00 


o 




1^ 


o 


t3\ 




m 




<N 




1^ 




00 






1^ 






a\ 










m 


IT) 






IT) 










<ri 








00 








td 


td 


td 


td 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 



o 



1^ 
d 



(N 00 



^^<ri<rim(N(Nooooi^^<ri 



o 



y3 



00 00 tX) 



1^ 

ON 



tON 

o 

t3N 



in M 00 ON M 
— ; m tN <N <N 

d 00 t3N ON d 



00 



o 



ON O 



ON 

00 d 



1^ -H 

ON d 



oo d 



IN <( 

a 8 

H c/1 



s ^ 

^ § 

o o 

o o 

o o 

o o 

o o 

o o 

•z. -z. 



Oh 

Q 

•z. 
u 



Ph 



CO oo 



U 

O 

u 

c/l 



§ i p 

tr 00 

^ M 

M m 



ION <J~i 



O W U P; 

J Ph Ph Oq O <f 
H .> Q X 5. 



^ NO 



s s a o 



U &o Pi 



hJ |jh U 00 



|Jh 

Ph S 

Pi 00 



S 



I X) ^ 

I r .J^. 



1^ 

NO 

o 



o 


<N 






1^ 




in 




00 




in 




m 


tON 


NO 


m 


o 








in 


00 


tON 


<N 






m 




NO 




rn 










tON 


NO 


in 


1^ 


1^ 










00 


oo 


1^ 


NO 


in 


00 




<N 


00 


m 


NO 


rn 




tON 


NO 


NO 


NO 


o 


NO 


NO 


rn 




tON 




00 


(N 


tON 


in 


'S- 


m 


in 


O 


O 


o 


o 


o 


o 


o 


o 




o 


o 








tq 








m 


(N 




<N 


!N 


(N 


m 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 



< 

NO 

u 

c/l 



NO 




in 




in 






in 


1^ 




C3N 


CNl 


in 


in 


in 




00 




NO 


1^ 


(N 


NO 


tON 


in 




in 


tON 


in 


00 




ON 


00 




00 


00 








m 








o 


!N 


00 


1^ 


NO 




00 






-a- 


in 


ON 


(N 


o 


00 






1^ 






in 


rM 


in 


(N 


in 


m 


m 


O 




in 




tON 


(N 


NO 


1^ 


o 


o 


1^ 


00 


NO 


m 


5^ 


00 


00 


ON 






o 


(N 


o 


in 


o 


o 


in 


m 


!N 


m 


NO 


00 




00 






00 




O 


(N 


in 


^ 


(N 


^ 


m 


o 


o 


o 


o 


o 


o 


O 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


O 


O 


o 


O 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


O 


o 


o 


o 


o 


o 


c/l 


c/l 


c/l 


c/l 


c/l 


c/l 


c/l 


c/l 


c/l 


c/l 


c/l 


c/5 


c/5 


c/l 


c/5 


c/3 


c/l 


c/l 


c/l 


c/l 


c/l 


c/l 


c/l 


•z. 


•z 


•z 


•z 


:z 




•z 


•z 


•z 


•z 


:z 




•z 


:z 






•z 






•z 


•z 




•z 




M 


M 


M 






M 


M 


M 


M 


M 


M 


M 


M 


M 


M 


M 


M 


M 


M 


M 


M 





2285 



Molecular Vision 2013; 19:2274-2297 <http://www.molvis.org/molvis/vl9/2274> 



) 2013 Molecular Vision 



73 



b > 



CUD 



o 

s 



o 
ft 



B 
o 



ft 
o 

ft 



1^ 








m 


NO 


o 




(N 


00 




ON 


o 


d 


d 




(N 


in 


00 


oo 


NO 


<ri 


ON 


IT) 


m 


o 


NO 


d 


d 


d 


m 


IN 


NO 






p 










(N 


o 












ON 



t< ON 



NO 



NO 



Oh Q ^ 

H E Q 



o 
o 
o 
o 
o 



•z. -z. 



2; 

M 



u 
o 

i 

o 



"3 
'3 



o 

S3 

o 
o 

M 

.a 

t3 



ft 

X 



ao 

T3 
U 

•K 
o 
ft 



ft 



2286 



Molecular Vision 2013; 19:2274-2297 <http://www.molvis.org/molvis/vl9/2274> 



) 2013 Molecular Vision 




Figure 6. Heatmap of genes highly expressed in the RPE. Dark shading indicates low expression; light shading indicates high expression. 



based on its sum of distances to other arrays (Figure 2B,D). 
However, as the magnitude of this difference was marginal, 
we retained this sample in the data set. After CTRL3 and 
CTRL9 were removed, the p value for age between the two 
groups was 0.70. 

Differentially expressed genes: AltAnalyze mapped the 
exon array probesets to 42,187 unique Ensembl identifiers. 
To reduce false positive observations, we filtered this set by 
selecting only protein coding genes with mean log^ expres- 
sion greater than or equal to the median expression of all 
protein-coding genes (7.11; Figure 3), which resulted in 10,286 
genes. We found 75 genes differentially expressed in AMD 
compared to controls (at least 50% difference in expression; 
moderated t test raw p value <0.01; Table 2; Figure 4). None 



of the genes were significantly differentially expressed after 
FDR correction (adj. p value <0.1). We submitted the 52 
downregulated genes and the 23 upregulated genes to DAVID 
for annotation. Of genes upregulated in early AMD, three 
genes (FADSl, FADS2, PTPLA) were identified as implicated 
in biosynthesis of unsaturated fatty acids (fold enrichment 
= 86.68; Benjamini corrected p value = 0.0052). Of genes 
downregulated in early AMD, 28 terms were significantly 
enriched (fold enrichment from 2 to 128; Benjamini-corrected 
p value <0.01). These terms were associated with vision, 
sensory perception, and the plasma membrane. No transcripts 
flagged by AltAnalyze as putative splicing events appeared 
to be alternatively spliced upon close visual inspection of the 
probeset data in DomainGraph. 
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Comparison with previous studies: Genome-wide asso- 
ciation studies have identified several genes harboring 
polymorphisms associated with AMD. We compiled a list 
of these genes and examined their expression levels in our 
data set (Table 3; Figure 5). Some genes (e.g., ADAMTS9, 
APOE, ARMS2, B3GALT1, CCR3, CETP, CFHR2, CFHR3, 
COLlOAl, DDR], LIPC, RAD51B, TGFBRl, TLR3) did not 
meet our filtering criteria (i.e., mean expression in AMD 
or control groups > median expression of all protein coding 
genes). Of the genes previously associated with AMD, none 
were differentially expressed (raw p value <0.01). 

Analysis of genes associated with either RPE or vascular 
endothelium: We next evaluated genes with cell-type-specific 



expression in either the RPE or endothelium. For RPE genes, 
we took the top 35 genes identified as highly expressed in the 
RPE versus the retina and choroid in the recent manuscript 
by Booij and colleagues [33]. We added SERPINFl (PEDF), 
another marker of RPE, to this set [34]. Of these, 28 mapped 
to Ensembl IDs in our filtered set and were included for 
analysis. Thirteen RPE-expressed genes showed between 
12% and 72% increased expression (raw p value <0.1; Table 
4; Figure 6). 

As a proxy for choroidal endothelial cells, we compiled 
a list of 24 genes with known expression in the endothelium 
and compared the fold change between the AMD and control 
samples. Endothelium-associated genes were based on gene 
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lists identified as specific to the endothelium [35] or from 
a literature search. Seven of these genes showed decreased 
expression in early AMD between 20% and 189% (raw p 
value <0.1; Table 5; Figure 7), which suggested decreased 
numbers of vascular endothelial cells in the choroid and/ 
or dedifferentiation of extant endothelial cells. Thus, there 
appeared to be a trend toward increased expression of RPE- 
expressed genes and decreased expression of endothelium- 
expressed genes in eyes with early AMD. 

To calculate the enrichment of the RPE and endothelium 
gene sets, we performed GSEA [29]. For each gene set, GSEA 
calculates (a) an enrichment score (ES) based on the rank 
distribution of individual genes within the set from among 
all unique genes symbols in our filtered set (n=10,286), (b) a 
normalized enrichment score (NES) that accounts for the size 
of the set, (c) a nominal p value based on phenotype permuta- 
tion, e.g., AMD or control status of samples, and (d) an FDR q 
value to control for the gene set size and multiple hypothesis 
testing. GSEA suggests a greater trend for the overall enrich- 
ment of the endothelium-associated gene set in AMD samples 
(ES = -0.78; NES = -1.39; nominal p value=0.133; FDR q 
value = 0.08) than the enrichment of RPE-specific genes in 
the control samples (ES = 0.81; NES = 1.3; nom. P value = 
0.167; FDR q value = 0.191). Since an FDR q value of less 
than 0.25 is considered significant for GSEA [29], these data 
suggest a significant decrease in endothelial cell transcripts. 

Analysis of CF// risk genotypes: Last, we reanalyzed the data 
set by stratifying based on genotypes at the rsl061170 SNP 
in CFH {W high-risk, YH/HH samples; five low-risk, YY 
samples), independent of AMD affection status. The same 
filtering criteria were applied as in the previous analysis, 
retaining 10,288 protein-coding genes. Thirty-five genes 
were identified as differentially expressed (at least 50% 
difference in expression; moderated t test raw p value <0.01; 
Table 6). None of the genes were significantly differentially 
expressed after FDR correction (p value <0.1). 

DAVID analysis of the 12 genes with increased expres- 
sion in the high-risk genotype revealed no significantly 
enriched terms (with Benjamini corrected p value <0.01). 
In the DAVID analysis of the 23 genes with lower expres- 
sion in the high-risk genotype, DAVID identified eight 
terms with significant enrichment (fold enrichment from 3 
to 8; Benjamini corrected p value <0.01). These terms were 
associated with extracellular secretion, signal peptides, and 
disulfide bonds. 
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DISCUSSION 

AMD is a complex disease that shows altered function and 
viability of photoreceptor cells, RPE, and choriocapillaris 
endothelial cells. Although our understanding of the genetics 
of AMD has progressed in the last decade, many basic ques- 
tions about the pathogenesis of this disease remain. In the 
current study, we evaluated gene expression at the exon level 
in AMD and control eyes and found a loss of transcripts 
expressed by choroidal endothelial cells. 

Limitations of the current study include small sample 
size (due to stringent inclusion criteria) and the lack of 
strong signatures of differential expression between the case 
and control samples, obviating validation of single genes. 
Increased sample size may reveal some genes with a small 
difference between the cases and controls, which remain 
significant after multiple hypothesis testing. However, this 
is unlikely. Previous analysis of AMD-affected tissues with 
microarray technologies required the use of non-standard 
methods to identify altered gene expression (expression corre- 
lation between two data sets with p<0.1 and 25% fold change) 
[15] or the use of sensitive clustering algorithms to identify 
patterns across many AMD grades simultaneously [14]. At 
the tissue level in our data set, the apparent trends among the 
RPE and endothelial genes (overall direction of fold change 
and p<0.1) are supported by the GSEA results. These gene 
set-level changes are small in our samples. 

Cross-contamination by photoreceptor-specific transcripts: 
We observed an elevation of select photoreceptor-specific 
genes (i.e., RHO, PDC) across all control samples compared 
to the AMD-affected samples, particularly noticeable in 
CTRL4 and CTRL5. High expression of retinal transcripts 
is commonly reported in gene expression studies of the RPE 
and choroid [14,36,37], and various strategies have been taken 
to deal with cross-contamination (e.g., laser capture microdis- 
section [33], flagging genes with high expression in one tissue 
that appear in another tissue [14]). Cross-contamination is not 
surprising, as photoreceptor cell outer segments are partially 
interdigitated and ensheathed by the apical microvilli of the 
RPE and supported by the interphotoreceptor matrix [38]. 
Although most photoreceptor cell transcripts are expected 
to be within the outer nuclear layer and inner segments. 
Van Soest and colleagues suggest that these transcripts 
may be present at the interface between the photoreceptor 
outer segments and apical RPE [37]. The photoreceptor cell- 
specific gene expression we observed could be an artifact 
of mechanical separation or a stochastic individual variation 
in adhesion, and this random event might have occurred to 
a higher degree in the control eyes. However, we hypoth- 
esize that the pattern of elevated retinal contamination in 
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the controls, which was absent from AMD cases, may be 
due to decreased adhesion of the retina to the RPE in the 
AMD samples compared to the controls. Although adhesion 
between the neural retina and the RPE has been examined in 
primates [39], to our knowledge this has not been investigated 
in the context of early AMD. Alternatively, if neural retina 
components are uniformly present in all samples, then the 
decreased expression of retina-specific genes in the AMD 
samples could indicate that expression within photoreceptors 
is decreased even at the earliest stages of AMD. 

Intriguingly, genes with lower expression in AMD than 
in controls included KDR and FLTl, known VEGF receptors. 
Hierarchical clustering reveals that these genes do not segre- 
gate with the neural retina-specific genes and that expression 
of this sub-cluster of genes is consistent across controls and is 
not affected by CTRL4 and CTRL5 (Figure 1). 

Loss of endothelial-specific gene expression: Our group 
has previously shown that the density of microvasculature 
decreases as the volume of sub-RPE drusen increases [9]. 
Similarly, studies of blood flow in eyes with drusen show 
susceptibility of the choroidal vasculature to degenerative 
changes in AMD [11]. In elegant whole mount studies of eyes 
with advanced AMD, McLeod and colleagues found a linear 
relationship between the RPE and choriocapillaris, with the 
loss of either layer affecting the integrity of the other [10]. 
From these studies, with a sample size similar to that in the 
current report, the authors concluded that RPE loss precedes 
endothelium loss in GA, whereas endothelium loss precedes 
RPE loss in CNV [10]. Our data indicate that gene expression 
of endothelium-related genes decreases in early AMD before 
any atrophy or neovascular change develops. 

Our results are also consistent with state-of-the-art 
proteomics studies of Bruch's membrane-choroid prepara- 
tions, in which the authors found the endothelial cell proteins 
von Willebrand factor (VWF) and carbonic anhydrase 4 
(CA4) are reduced in early and mid-stage AMD and advanced 
dry AMD, respectively [40]. The same authors showed no 
loss of RPE proteins RPE65 and CRALBP (RLBPl). In 
fact, the levels of these proteins were elevated in Bruch's 
membrane, consistent with our findings at the mRNA level. 
Taken together, these studies suggest loss or dedifferentiation 
of choroidal endothelial cells before loss of the RPE in eyes 
with AMD. In addition, expression of the choroidal endothe- 
lial genes ICAMl, SELE, and PLVAP was decreased in some 
AMD classes in the report by Newman et al. (supplemental 
data) [14]. 

Decreased expression of ADAMTS9 in samples with high 
AMD risk CFH genotype: Eyes with high-risk CFH geno- 
types were also assessed for gene expression changes. We 
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previously found increased membrane attack complex 
formation in eyes homozygous for the high-risk allele [41], 
and complement activation may affect gene expression in 
nearby cells. Although not significant after multiple hypoth- 
esis testing correction, expression of ADAMTS9 transcripts 
was lower in samples with high-risk CFH genotypes than in 
samples with the low-risk genotype. This gene is a member 
of the ADAMTS (a disintegrin-like and metalloprotease 
domain with thrombospondin type I motifs) protein family 
with the ability to cleave aggrecan and versican [42]. Recently 
a SNP (rs6795735) located 32.5 kb upstream of the ADAMTS9 
transcription start site was associated with increased risk of 
AMD in a meta-analysis of genome-wide association (GWA) 
data [43]. ADAMTS9 is expressed in ARPE-19 cells [44] and 
several types of microvascular endothelial cells [42], although 
to our knowledge localization of ADAMTS9 has not been 
examined in human RPE and choroid samples. In microvas- 
cular endothelial cells, ADAMTS9 suppresses angiogenesis, 
albeit not via sequestration of VEGF165 [42]. Further experi- 
mentation is necessary to elucidate the relationship between 
CFH activity and ADAMTS9 function in early AMD. 

In conclusion, we performed microarray analysis of 
human donor maculas and found early loss of choriocapillaris 
endothelial cell markers in early AMD, with preservation of 
RPE cell transcripts. These results have potential impor- 
tance for therapy. Recently, the replacement of RPE cells 
in AMD has been contemplated as a treatment for AMD. 
These studies, as well as proteomic [40], anatomic [9,10], 
and clinical [11] studies, strongly suggest caution in this type 
of approach, since transplanting healthy RPE into a macula 
with a degenerated choriocapillaris may be fruitless. In addi- 
tion to replacing photoreceptor cells and RPE, strategies for 
replacing lost choriocapillaris are necessary to fully restore 
function in AMD. 
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